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Abstract. We discuss the results of the MEGA microlensing campaign towards M31. Our analysis is based on 
an analytical evaluation of the microlensing rate, taking into account the observational efficiency as given by 
the MEGA collaboration. In particular, we study the spatial and time duration distributions of the microlensing 
events for several mass distribution models of the M31 bulge. We find that only for extreme models of the M31 
luminous components it is possible to reconcile the total observed MEGA events with the expected self-lensing 
contribution. Nevertheless, the expected spatial distribution of self-lensing events is more concentrated and hardly 
in agreement with the observed distribution. We find it thus difficult to explain all events as being due to self- 
lensing alone. On the other hand, the small number of events does not yet allow to draw firm conclusions on the 
halo dark matter fraction in form of MACHOs. 
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1. Introduction 



Since the proposal of IPaczvnskil l|l98rJ) gravitational mi- 
crolensing has probed to be an efficient tool for the study 
of the MACHO contribution to the dark matter galac- 
tic halos. The first line of sight to be explored with this 
purpose has been that towards the Magellanic Clouds 
jAlcock et all Il993t lAubourg et all 11993 lUdalski et all 
Il993l) As fir st discu s sed by ICrottsf |l 992(1 MBaillon et all 
jl993[) and IJetzeJ ()l994f) o bservations toward s M31 
have also been u ndertaken ijCrotts k, Tomanevl Il996l 
lAnsari et ailll997|) . 

The interpretation of the results obtained so far re- 
main, however, debated and controversial. Along the line 
of sight towards th e LMC the MACHO collaboration 
l|Alcock et alll2000|) reported the signal of a halo fraction 
of about 20% in form of MACHOs with mass ~ 0.5 M Q , 
while the latest results of the ER OS collaboration to- 
wards both the SMC and the LMC llAfonso et all 120031 
iTisserandl 120051 iTisserand et"afl |2006) are even compat- 
ible with a no MACHO hypothesis. 

The case towards M31 is complicated by the further 
degeneration in the lensing parameter space due to the 
fact that sources at baselin e are unresolved, a case re- 
ferred to as " pixel-l ensing" ijCrottsl Il992t [Baillon e t all 
Il998t iGouldl Il99f^ . Still, a handful of microlensing 
event s have been obs er ved in the meantime l|Ansari et ah , 
Il999t lAuriere et all 120011 ICalchi Novati et all l20~ , 



20031 IPaulin- 


Henriksson et all 120031 iRiffeser et all 12002 


de Jong et al. 


. 12004 iBelokurov et all 120051 IJoshi et al. 



have been rep orted. The results of the POINT- AGAPE 
collaboration l|Calchi Novati et all l2005|) are compatible 
with the ones of the MACHO group, by putting a lower 
limit on the halo fraction in form of MACHOs of ~ 20% 
for objects in the mass ra nge 0.5 — 1 On th e contrary 
the MEGA collaboration l|de Jong et all 12006^1 finds that 
their results, although not conclusive, are in agreement 
with a no MACHO hypothesis. Although the issues in- 
volved in the microlensing observations towards the LMC 
or the M31 are indeed rather different, the results on the 
halo fraction in form of MACHOs depend crucially on the 
prediction of the expected signal due to known luminous 
populations, this being dominated by the "self-lensing" 
signal where both source and lens belong to same star 
population residing respectively either in the LMC or in 
M31. This problem is indeed the main aspect we want to 
discuss in this paper. 

The issue of the expected microlensing si gnal towards 
M31 has been discussed in a few works (e.f 
20011 iBaltz et al.ll2003t IRiffeser et all 12006 



Kerins et al 



Kerins et al 



2006). In this respect the modeling of the M31 luminous 
components is a crucial aspect to be dealt with in order 
to get me aningful results. We ha ve first considered these 
aspects inlDe Paolis et ahl l|2005|) . then, more recently, in 
llngrosso et al.l l)2006[) we have developed a Monte Carlo 
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simulation which we have used to investigate the nature 
and location of the microlensing candidates events towards 
M31 a s reported in a first paper by the MEGA collabo- 
ration l|de Jong et all l2004j) . In the present work our aim 
is to further explore t hese issues taking into account the 
latest MEGA results l|de Jong et all 12006^ . In particu- 
lar, we now go through a full characterisation of the ex- 
pected signal, including the predicted number of events, 
that we then compare with the observational results. Our 
aim is to explore, in particular, the question whether the 
expected self-lensing signal due to stars belonging either 
to the bulge or the disc of M31 is able, as claimed by 
Ide Jong et al] (|2006|h to fully explain their results. 

The plan of the paper is as follows. In Sect. 2 we de- 
scribe the microlensing rate, our main tool of investiga- 
tion, and present its predictions. In Sect. 3 we critically 
discuss the models used to describe the different galac- 
tic components involved. In Sect. 4 we discuss our main 
results and give some concluding remarks. 



2. Event rate prediction 

In evaluating 1 the expected event number along a fixed 
line of sight we take into account the existence of two 
source populations (stars in the M31 bulge and disk) with 
number density n s (D os ,M) and of five lens populations 
(stars in the M31 bulge, stars in the M31 and MW disks, 
MACHOs in M31 and MW halos) with density n x (D ol ^i). 
Here D os (D Q i) is the source (lens) distance from the ob- 
server, M is the source magnitude and /i is the lens mass 
in solar units. 

We assume, as usual, that the mass distribution of 
the lenses is independent on their position in M31 or in 
the Galaxy (factorization hypothesis). So, the lens number 
density (pe r unit of volume a nd mass) ni(D i,ii) can be 
written as ijjetzer et all 2002) 



ni{D oh fi) 



PljDol] 

Mo) 



■00 (m) 



(1) 



where pi(0) is the local mass density of the considered 
lens population in the Galaxy or the central density in 
M31, ipoif 1 ) the corresponding lens number density per 
unit mass and the normalization is given by 



Mup 



00 (A 4 ) t 1 d V 



MO) 

Mr, 



(2) 



Here /i m i n and /z up are the lower and the upper limits for 
the lens masses (see Subsection 3.2). 

Likewise, assuming that the magnitude distribution of 
the sources is independent on their position in M31, the 
source number density (per unit of volume and magni- 
tude) n s (D os , M) can be written as 



(3) 



1 Here we follow with some modificat ions the derivation in 
our previous paper Irngrosso et alJ i2006T) . 



Table 1. For the 14 MEGA events we give position, mag- 
nitude at maximum Ar and full-width half-maximum du- 
ration ii/2- The coordinate system we adopt has origin in 
the M31 center and the X axis oriented along the M31 
disk major axis. 



MEGA 


X 


Y 


Ar 


tl/2 




arcmin 


arcmin 


mag 


day 


1 


-4.367 


-2.814 


21.8 + 0.4 


5.4 + 0.7 


2 


-4.478 


-3.065 


21.51 + 0.06 


4.2 + 0.7 


3 


-7.379 


-1.659 


21.6 + 0.1 


2.3 ± 2.9 


7 (N2) 


-21.164 


-6.248 


19.37 + 0.02 


17.8 + 0.4 


8 


-21.650 


+7.670 


22.3 + 0.2 


27.5 + 1.2 


9 


-33.833 


-2.251 


21.97 + 0.08 


2.3 + 0.4 


10 


-3.932 


-13.847 


22.2 + 0.1 


44.7 + 5.6 


11 (S4) 


+19.193 


-11.833 


20.72 + 0.03 


2.3 + 0.3 


13 


+22.072 


-22.022 


23.3 + 0.1 


26.8 + 1.5 


14 


+19.349 


-29.560 


22.5 + 0.1 


25.4 + 0.4 


15 


-6.634 


-0.697 


21.63 + 0.08 


16.1 ± 1.1 


16 (Nl) 


-6.886 


+3.843 


21.16 + 0.06 


1.4 + 0.1 


17 


+21.214 


-5.161 


22.2 + 0.1 


10.1 ± 2.6 


18 


+6.995 


-13.533 


22.7 + 0.1 


33.4 + 2.3 



where £ s (0) is central luminosity density of the considered 
source population, <j) s (M) is the source number density per 
unit magnitude in the M31 center and the normalization 
now reads 



<j) s (M) L(M) dM =£ s (0) 



(4) 



Here M min and M up are the lower and the upper limits 
for the source magnitude (see Subsection 3.1), L(M) is 
the luminosity in a given photometric band 



L{M) = rjvega L Q 10- M / 2 - 5 



(5) 



VVega being the Vega luminosity (in solar units) in the 
considered band. 

We consider the volume element of the microlensing 
tube to be d 3 x — (v t ± ■ n)REUthdadD \, Re being the 
Einstein radius, v t ± the relative tranverse velocity be- 
tween the lens and the microlensing tube with distribu- 
tion function f(v r ±), n the inner normal to the microlens- 
ing tube and a the angle between n and A± (see eq. |HJ- 
Assuming perfect observational sensitivity to microlens- 
ing, the differential event rate dN cv /dfl (in units of event 
sr _1 ) for microlensing by compact objects with impact pa- 
rameter below a cert ain threshold Uth, during the time in- 
terval dt, is given by ijGriestl Il99ll be Ruiula et all Il99lh 



dN cv 

dn 



D 2 os u th R E da v r± f(v I _ L )d 2 v I± cos9 (6) 
ni(D Q i,fi) n s (D os ,M) d/i dM dD os dD \ dt , 



where 9 E (— 7r/2,7r/2) is the angle between n and v r ±. 
We assume that the velocity distributions of lenses and 
sources are isotropic around their streaming velocities (if 
present) due to the rotation of the considered popula- 
tion with respect to the M31 or MW center (we neglect 
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any transverse drift velocity of the M31 center with re- 
spect to the Galaxy). Accordingly, the lens (source) ve- 
locity is splitted into a random component - which fol- 
lows a Maxwellian distribution with one-dimensional ve- 
locity dispersion <j\ (<r s ) - and a streaming component, 
namely V\ = t>K ran + v Lmt and v s = t> Siran + v s , mt . When 
the lens and source velocities are projected in the lens 
plane (transverse to the microlensing tube), the respec- 
tive random velocity distributions are again described by 
Maxwellian functions, with the same one-dimensional ve- 
locity dispersion o\ for lenses, and with (projected) disper- 
sion {D\/D s )a s for sources. Then, neglecting the stream- 
ing, the relative, projected, random velocity t>i s ±,ran = 
■wii.ran ~ w s_L,ran of lenses and sources is a maxwellian dis- 
tribution / (uisx ran) with combined width 



<7 S 1 = J<J 2 + (D ol /D os ) a 2 



(7) 



We now include all streaming motions in the vector A± 
defined as the difference between the projected, streaming 
velocities of lenses, sources and observer, namely 

D ol \ ^f D ol\ ,„x 

1 - £j- ) v Q±,rot +[ — ) V ^rot ~ WlX.rat -(8) 



D n 



The resulting distribution function f(v T ±) of the rela- 
tive, transverse velocity between the lenses and the mi- 
crolensing tube is now given by the Maxwellian function 
/(i>ls-L,ran) shifted by the vector A±, that we write in po- 
lar coordinates on the lens plane as 



f(v I± )d 2 v r± = 



1 



t! r _i_ dv t ± d9 



(9) 



Taking a to be the angle between A± and the normal n 
to the microlensing tube, it results that ip = a + 9, where 
<p is the angle between v r ± and A± . 

We recall that in the pixel lensing regime the effective 
radius of the microlensing tube is a function of the source 
star magnitude, namely Uth = u t h(M). Moreover in the 
following we evaluate the differential rate taking into ac- 
count an efficiency function that depends on the impact 
parameter, e = e(uth)- Therefore, we are going to replace 
in eq. JBJ dN cv by J dN cv /duth x e(it t h)duth, with upper 
limit ut(M). 

Eventually, after integration on the angular variables 6 
and a, one obtains the expected event number rate (events 
sr _1 ) during the observation time T OD s 



dN cv r- /4GM Q 

= Jobs 4V2<7si' 



u T (M) 



du 



ih 



M up rfJ-up poo 

4> s (M)dM / dn^Mv) / D 2 0S dD 

A/min Jflmin JO 



(10) 



D 



dD. 



oH 



I D os { A (0) J \ £ B (0) 

:+ ^I (2f3z) e(t 1/2) A/) dz . 



where z = v I j_/(^/2a s \), (3 = \ A±]/{^2a s \) and Iq{2(3z) is 
the zero-order modified Bessel function 2 . 

In the previous equation we explicitely take into 
account an experimental event detection efficiency 
e(ti/ 2 ,A/), given as a function of the full-width half- 
maximum event duration ti/2 



h/2 = t E f(a) , 

/(a) = 2V§(-^ 
\Va 2 



a + 1 



4a i/a 2 + 2a 



1/2 



(11) 



and of the maximum flux difference during a microlensing 
event 



A/ = f (A n 



1) 



(12) 



Here is the Einstein time, A max = A max (uth) the am- 
plification at maximum and /o the unlensed source flux. 

It is well known that self-lensing and dark-lensing 
events may have different time durations, depending on 
the MACHO mass value. On the other hand, in pixel- 
lensing observations experimental results are usually given 
in terms of the ti/2 time scale. Thus, it is important to 
evaluate the expected event rate as a function of ti/2- 

From eq. (|ll|l and the relation t E = Re/vi± it follows 



n/2 



R E f(a) 
zV2ct s i 



(13) 



and it is straightforward to derive the differential event 
rate 

ut(M) 



^cm V (*Va) = T ° bs 8cr si / du th 



M up /"Mup roc 

4> s (M)dM / d[iMp) / D 2 0S dD 

Afmto Jfhmlx, JO 



(14) 



where z is now given in terms of and A max through 
eq. jnj). 

The model parameters that need to be specified are 
the luminosity <fi s (M) and mass ipo(l^) functions, the stel- 
lar mass distributions in M31 and MW, the mass-to- 
luminosity ratios for the stellar populations in M31, the 
velocity dispersion cr s and o\ for the source and lens popu- 
lations. Further model parameters derive from the consid- 
eration of the existence of dark matter in both M31 and 
MW halos. 



2 By comparing eq. 111011 with eqs. (11) and (12) in 
llngrosso et alJ (I2006T) one can see that the composition of the 
two maxwellian (projected) velocity distributions for lenses 
and sources permits now to evaluate analytically the two- 
dimensional integration on the source velocity in eq. (12). 
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Table 2. The M31 disk and bulge models. Relevant parameters for WeCapp l|Riffeser et alll2006l) and our reference 
model. Values of mass density, distance, mass and velocity are given in units of M pc~ 3 , kpc, 10 10 M Q and km s -1 , 
respectively. For the reference model, the symbol — means that the corresponding value as in the WeCapp m odel is 
used. In the last row we give some relevant parameter values for the models MEGA A in lde Jong et al.l 1 2006F) . 





bulge 


disk 


model 


p(0) a M ext B (M/L R ) a 


p(0) h H M ext fl (M/L R ) a(2h) 


WeCapp 
reference 
MEGA A 


3.97 x 10 4 4.00 0.36 2.96 140 
4.53 x 10 4 2.62 x 10" 3 3.85 

4.4 3.6 


0.20 6.4 0.3 3.09 0.68 0.88 40 
1 5.5 2.4 



3. Models 

3.1. Source luminosity function 

In pixel-lensing experiments only bright and sufficiently 
magnified sources can give rise to detect able microlens- 
ing ev ents. Monte Carlo simulations (e.g. Ilngrosso et alJ 
( 200(|)) allow to determine the useful range of source mag- 
nitude M m j n ~ — 6 and Af max ~ 3, and the threshold value 
for the impact parameter ut(M). 

As regards the source luminosity function 4> S (M), in 
the lack of precise information about the luminosity func- 
tions in M31, we adopt the luminosity function derived 
for local stars in the Galaxy and assume that it also holds 
for M3 1, irrespectively on the pos ition. In particular, fol- 
lowing |Mamo^^n^^onei3 l|l982(l . the stellar luminosity 
function in the magnitude range —6 < M < 16 is given 

by 

MM) = H [1 + 1Q _ (a _ mM _ M , )]1/s , (15) 

where, in the i?-band (the observational band of the 
MEGA collaboration) M* = 1.4, a ~ 0.74, (3 = 0.045 
and 5 = 1/3. The constant H in eq. (|15|) is determined 
via the normalization condition in eq. (@J , namely 

j <f> s {M) L(M) dM = M0) f|£) , (16) 

where (M/Lji) is the mass-to-luminosity ratio for the 
source star population in the i?-band. Note that the nor- 
malization for the source density distribution, eq. I|20|) . 
implies that the event rate does not depend on (M/Lr). 



mean mass for lenses in the bulges and disks are (mj) ~ 
0.41 M & and (mA ~ 0.51 M Q , respectively. 

We also consider steeper mass function as proposed by 
Zoc cali et al. 1 1I2OOOI) and wc find that our estimate of the 
self-lensing event number turns out to be rather insensitive 
to the mass function choice. 

For the lens mass in the M31 and MW halos we assume 
the (S-function approximation 

M,) = (18) 
Mil 

and take a MACHO mass, in solar units, /ih = 
10" 1 , 0.5, 1. 

3.3. Mass distributions in M31 and MW 

The visible mass distributions for the M31 bulge and disk 
ar e derived by fi tting the observed brightness profiles given 
by iKentl l)l989|) and by further assuming mass- to-light ra- 
tios for bulge and disk stellar populations. Moreover, the 
consideration of the M31 rotation curve data allows us to 
derive the distribution of the dark matter in the M31 halo. 

Here, we use a coordinate system (x,y, z) centered in 
M31, with x axis along the disk major axis. We also as- 
sume that the disk is inclined by the angle i = 77° and 
that the disk azimuthal angle relative to the near minor 
axis is 4> = 38.6°. The position angle of the bulge is 50°. 

We neglect the MW disk since we have verified that 
the expected number of events due to lenses belonging to 
this mass component is only about 1% of the total number 
of M31 self-lensing events. 



3.2. Lens mass function 

As far as the lens mass function is concerned, for lenses 
belonging to the bulge and disk star populations, the 
lens mass is assu med to follow a broken power law 
JOould et alllTflflj 

Mm) = #1 M~ 56 for Mmin < fi < 0.59 

= K 2 fi" 2 - 20 for 0.59 < fi < /i U p (17) 

where the lower limit /i m i n = 0.1 and the upper limit /x up 
is 1 for M31 bulge stars and 1.7 for M31 and MW disk 
stars. The constants K\ and K2 are fixed according to the 
normalization condition given by eq. (J5J). The resulting 



3.3.1. M31 bulge 

The M31 bulge model is derived from Tab. I in lKentl <|l989f) 
containing the bulge 3-d brightness density in the Gunn 
r-band and the ellipticity e(a) as a function of the major- 
axis distance a to the M31 center. 

We fit the 3-d brightness profile with a single de 
Vaucouleurs a 1 / 4 law (reference model) 

Jr(o) = J,(0) io-°- 4 ( 7 - 598 - 1/4 ) a > a min , (19) 

with central 3-d brightness density j r (0) = 9.57 x 10~ 7 
Lq arcsec~ 3 (shifting to magnitudes, eq. (|19(l may be writ- 
ten in the form m r (a) = 15. 048 + 7. 598a 1 / 4 mag arcsec" 3 ). 
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This model accurately fits Kent data for a mm ~ 1 arcmin, 
namely in the region usually explored by pixel lensing ob- 
servations. 



1 4 r 




0.01 



0.10 



.00 10.00 
a (arcmin) 



100.00 



1000.00 



Fig. 1. The projected 2-d brightness profile (solid line) is 
shown for the reference model in comparison with Kent 
data (crosses). Dotted and dashed lines give the bulge 
and disk contribution, respectively. The dotted dashed line 
shows the bulge contribution for the boxy model. 



has the same behaviour as the brightness profile and cen- 
tral mass density given by 



P(0) 



M 



^Q-O.4[15.O48-(r-fi)-ext K -M R-d m od] (20) 



where (M/Lr) is the mass-to-light ratio in the i?-band, 
(r — R) the color of the bulge stellar population, Mqr = 
4.42 the absolute brightness of the Sun in the i?-band, 
extfl, the extinction in the same filter and distance mod- 
ulus rf mod = 24.43 (for an M31 distance of 770 kpc). By 
using the values (M/L p) = 2.96, (r-R) = 0.59 and 
cxt fi = 0.36 quoted bv iRiffeser et alJ £2006), we obtain 
p(0) = 4.53 x 10 4 Mq/pc 3 , corresponding to a total bulge 
mass Mbuigc — 3.85 x_10 10 M© , in agreement with the 
value given hv lKentl ill 9891). 

Note that the observed 2-d brightness profile is also 
compatibl e with more concent rated mass distributions for 
the bulge ((Beaton et alllioOfil) . For instance, we have tried 
a (boxy) model with 99% of the total mass inside 17.86 
arcmin (4 kpc). The mass density is now given by 

p(a) = 4.40 x 10 4 i -°- 4 ( 7 - 598a ° 24 > a < 17.86' 

= 1.81 X l 39 10 -O.4(7.598a o ») fl > 1? gg, _ (2J) 

In Fig. ^the projected 2-d brightness profile (in units 
of mag arcsec -2 ) is shown for both models together with 
Kent data. In deriving these profiles, we assumed that the 
bulge isophote are triaxial ellipsoids with semi-major axes 



300 



250 - 



200 - 



150 - 



100 - 



50 




50 



100 
a (arcmin) 



150 



Fig. 2. The full M31 rotation curve (solid line) is shown 
in comparison with data points d erived from HI mea- 
surements of Brinks and Burtonl l) 19841) (crosses) and 



ICarignan et al 



(2006) (diamonds). Dotted, dashed and 
dot-dashed lines give the bulge, disk and halo contribu- 
tion, respectively. 



From the 3-d brightness density profile in eq. Ijl9|l one 
can derive the corresponding mass density profile, which 



a 2 (e) 



V 



{1-ef 



(22) 



and ellipticity varying on the semi-major axis according 
to the Kent data 3 . From Fig. ^ on e can see that beyond 
0.03 arcmin both reference and boxy models accurately 
reproduce Kent data. The only difference is the behaviour 
of the bulge contribution at large distance where in any 
case the disk contribution is dominant. 

For comparison we also discuss the results obtained by 
using the bulge model ado pted by the WeCapp collabora- 
tion ( Riff eser et all Ef 



3.3.2. M31 Disk 

As in lKerins et all l|200l[) . the disk 3-d brightness density 
in the r-band is modeled by the law 



j r (x,y,z) = j r (0) eM-V* 2 +V 2 /h) sech 2 (z/if) , (24) 

and a best fit procedure to the Kent data (for a > 6 

arcmin) allows to obtain the central brightness density 
j r (0) = 4.2 x 10~ 13 Lq arcsec -3 (corresponding to a cen- 
tral magnitude m r (0) = 20.5), the radial scale length 

3 The existing relation between e and a may be approximated 
by iRiffeser et al.Ll200d) 



1 - e(o) 



0.254- 



+ 1.11 



(23) 



6 



Ingrosso et al.: Microlensing events towards M31 



h = 27.95 arcmin and the vertical scale length H = 1.34 
arcmin (corresponding to h — 6.4 kpc and H — 0.3 kpc, 
respectively) . 

As for the bulge, the corresponding disk mass den- 
sity profile follows the same behaviour as the bright- 
ness profile. Accordingly, the disk central mass density 
is derived by assuming the following parameter values 
(M/ Lr) = 0.88, (r - R) = 0.54 and ext R = 0.68 for the 
disk l)Riffeser et all l2006j) . implying p(0) = 0.2 M Q pc" 3 
and a total disk mass M ~ 3.09 x 10 10 M Q . The 2-d disk 
brightness profile is also shown in Fig. (dashed line). 

3.3.3. M31 and MW halos 

Both M31 and MW halo mass distributions are modeled 
as isothermal spheres 




p(r) = 



Po 



(25) 



For M31 a fit to the M31 rotational curve by using 
the three (bulge, disk and halo) component model al- 
lows to get the best fit parameter values rn = 2 kpc 
and p(0) = 0.23 pc~ 3 (see also iKerins et all l(200l|) 
and iRiffeser et all l|2006l) 'l. The overall M31 rotational 
curve and the contributions of the three components is 
shown in Fig. [2 In comparison with the recent determi- 
nation of the mass distribution in M31 ICarignan et alJ 
(2006), we find that at R = 35 kpc the dark mat- 
ter mass is — 3.7 x 10 11 M & and the stellar mass 
M vis = 6.6 x 10 10 M Q . This translates in a total dynami- 
cal mass of ~ 4.4 x 10 11 Mq and in a rotational velocity of 
233 km s _1 at R = 35 kpc, in agreement with the recent 
observations. The M31 halo is truncated at R — 150 kpc. 

For the MW we use a core radius a ~ 5.6 kpc and a 
local (i?o = 8.5 kpc) dark matter density p(Rq) — 1.09 x 
10 7 Mq kpc -3 . The corresponding asymptotic rotational 
velocity is u ro t — 220 km s _1 . The MW halo is truncated 
at R ~ 100 kpc. 



Table 3. The MEGA event detection efficiency 
e(ti/2, A/) is given as a function of 1/A/ (first row) for 
different values of tin (first column) in days. The nu- 
merical values are derived from Fig. 12 in Ide Jong et all 
J2£ 



Fig. 3. The map dN cv /dfl of the expected (total) event 
rate towards M31 is shown, assuming the reference model, 
a MACHO mass value p h = 0.5 and a MACHO halo 
dark matter fraction /h = 0.2. Here and in the follow- 
ing figures and tables we adopt the observational parame- 
ters of the MEGA collaboration. Accordingly, we consider 
T Q bs = 2 yr and we account for the detection efficiency 
e(ti/2)A £) and maximum im pact parameter ut(M) as 
given bv lde Jong et al.1 ((2006). From the outer to the in- 
ner M31 
5x 10 
arcmin 



region, contour levels correspond to the values 



3 , 1 x 10~ 2 , 2x 10~ 2 
2 respectively. 
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1 x 10 event 
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5 
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0.01 




10 


0.29 


0.30 


0.30 


0.19 
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0.01 
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0.25 


0.24 


0.20 
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0.05 


0.02 
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0.14 


0.16 


0.22 


0.21 


0.19 


0.15 


0.12 


0.08 



Fig. 4. The same as in Fig. [3] but for the dark-to-total 
event number ratio. From the inner to outer region, con- 
tour levels correspond to the values 0.4 , 0.5 , 0.6, 0.7 and 
0.8, respectively. 



3.4. Velocity dispersions 

The random velocities of stars and MACHOs are assumed 
to follow Maxwellian distributions, with one-dimensional 
velocity dispersion a = 140 and 166 km s _1 for the 
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Table 4. The integrated number of expected events inside 
each iso-rate contour of Fig.^is given for self-lcnsing and 
dark-lensing, assuming the reference model with fi^ = 0.5 
and /h = 0.2. In the last column we show the correspond- 
ing number of events detected by the MEGA collabora- 
tion. 



Events inside 










the 8 MEGA 










fields 












iso-rate contour 


self 


dark 


MEGA 




event arcmin -2 










1 x 10" 1 


3.80 


0.90 


1 




3 x 1(T 2 


6.49 


2.61 


4 


reference 


2 x 1(T 2 


7.45 


3.79 


5 




1 x 10" 2 


8.60 


6.77 


6 




5 x 10" 3 


9.20 


9.68 


12 




overall 


9.68 


11.76 


14 



M31 bulge and MACHOs, and a = 156 km s" 1 for 
the MACHOs in the MW halo. Moreover, following 
l|Widrow and Dubinskil l2005h . the M31 disk stars are as- 
sumed to have one dimensional dispersion velocity de- 
creasing towards the outer part from the central value 
cr(r = 0) ~ 110 km s" 1 to a{r = 30 kpc) ~ 5 km s" 1 . 
In addition, a rigid rotational velocity of 4 km s" 1 has 
been taken into accou nt for the M31 bulge l|Kerins et all 
1200 lUAn etall l2004h . For the M31 disk component the 
full rotational velocity as shown in Fig. (solid line) is 
also considered. 

4. Results and concluding remarks 

The main purpose of the present analysis is to compare 
the predictions of our model with the o bservational re- 
sults obtained by the MEGA collaboration l|de Jong et all 
200rj|. Therefore, to evaluate the microlensing rate we re- 
produce the MEGA observational set up and we make use 
of the event detection efficiency e(£]/2; A/), as a function 
of the time duration and amplification at maximum and 
th e maximum impact p arameter ut(M) values as given 
bv lde Jong et all lj200m . In Tab.El we give typical detec- 
tion effi ciency values derived from Fig. 12 in lde Jong et alJ 
(2006). In order to take into account the spatial varia- 
tion of the detection efficiency we use two different evalu- 
ations of e at distances smaller and larger than 11 arcmin 
from the M31 center 4 . It results that on average e is re- 
spectively smaller and larger by about 30% of the values 
quoted in Tab. [3] 

In the following tables and figures, we assume 
for both M31 and MW halos a MACHO halo dark 
matter fraction /h = 0.2 as suggested by mi- 
crolensing ob s ervati ons towards the Magellanic Clouds 
I Alcock et all l2000l) and pixel-lensin g observations to- 
wards M31 l|Calchi Novati et all l2005(K However, most of 
our results can be easily rescaled to other values of f\ x . In 

4 de Jong, private communication. 



Tab.|Hlwe consider different values for the MACHO mass: 
jUh = 0.1 , 0.5 and 1 (in solar units). Figures 3-6 and 
Tabs. EICI and |H1 are given for ^ = 0.5. 

Assuming the reference model for the M31 mass distri- 
bution, the spatial distribution of the expected events is 
shown in Figs.|3]and^ Here we give maps in the sky plane 
of the (total) event rate and dark-to-total event number 
ratio, respectively. In Tab. 0] we give our estimate of the 
integrated number of expected events inside each iso-rate 
contour of Fig. |2| Here and in the following we consider 
events inside the 8 fields selected by the MEGA collab- 
oration (as reported in Fig. 15 in Ide Jong et alJ (|2006) 
the innermost M31 region is excluded). From Fig. [3] and 
Tab. 01 one can see that dark-lensing gives an important 
contribution to pixel-lensing beyond the second (from the 
inner) iso-rate contour, namely beyond ~ 10 arcmin from 
the M31 center. 

The expected number of self-lensing events inside the 
8 MEGA fields is given in Tab. for different source and 
lens populations. Here with the symbols b, d and h we 
indicate sources and/or lenses in the M31 bulge, disk and 
halo, respectively. Capital symbol H is used to indicate 
lenses in the MW halo. In any case, the first (second) 
symbol refers to the source (lens). From Tab. 03 one can 
sec that for all the considered models (reference, boxy and 
WeCapp) the total number of self-lensing events is roughly 
the same (within 15%). 

As far as the reference and boxy models are concerned, 
we note an increase of bulge-bulge events to compensate 
a decrease of disk-bulge ones. This is expected to be due 
to the different concentration of bulge mass for the two 
distributions. We also note the increase of the disk-bulge 
events in the WeCapp model due to the more extended 
bulge mass distribution. 

Table 5. Number of self-lensing events expected given the 
set up of the MEGA campaign, for the different models 
discussed in the text. We consider different source and lens 
populations. 



Events inside 
the 8 MEGA 
fields 






bb bd db dd self 


reference 


4.25 1.17 3.30 0.96 9.68 


boxy 


5.14 1.10 2.76 0.95 9.95 


WeCapp 


4.98 1.34 4.08 0.96 11.37 



Assuming the reference model and /h = 0.2, in Tab. 
we give our estimate of the expected number of dark- 
lensing events for several MACHO mass value. We find 
that the total number of dark-lensing and self-lensing 
events turns out to be roughly the same. As regards the 
total (self+dark+background) number of expected events, 
~ 23 including ~ 1 event due to SN contamination (see 
next), it is consistent at 2tr confidence level with the 14 
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Table 6. For the reference model, the expected number 
of dark-lensing events is given for /ih =0.1, 0.5, 1 and 
/h = 0.2. 



Events inside 
the 8 MEGA 
fields 








Hh 


bh bH dh dH dark 


reference 


0.1 
0.5 
1 


2.55 1.04 8.81 3.10 14.49 
1.96 0.72 6.85 2.23 11.76 
1.68 0.58 5.80 1.82 9.88 



Table 7. The number of self-lensing and dark M31- 
lensing (due to MACHOs in the M31 halo) events is given 
for our two models (here labelled reference A and boxy A) 
assuming M b = 4.4 and M d = 5.5 (in units of 10 10 M©) 
and /ih = 0.5, /h = 0.2. To have the same luminosity for 
the M31 bulge and disk here we take (M/Lr)i, = 3.38 and 
(M/L R ) d = 1.56. We refer to models with M d = 5.5 and 
H = 1 kpc as maximal disk mode ls. In the last row w e 
report some results from Tab. 5 in lde Jong et alJ ((2006) , 
for the MEGA models in the case of high (MEGA A2) and 
low (MEGA Al) extinction and for a 20% M31 MACHO 
halo. 



Events inside 










the 8 MEGA 


self 


dark 


self 


dark 


fields 




M31 




M31 




H = 0.3 


H = 0.3 


H = 1 


H = 1 




(kpc) 


(kpc) 


(kpc) 


(kpc) 


reference A 


12.4 


8.8 


15.5 


8.6 


boxy A 


12.7 


8.5 


15.5 


8.7 


MEGA A2 






12.4 


5.7 


MEGA Al 






14.2 


6.2 



candidate MEGA events assumed to follow a Poisson dis- 
tribution. 

A comparison of our result s with the cor r espon ding 
values reported in Tab. 5 of Ide Jong et a l. ( 2006) for 
low and high internal extinction shows that there is a 
fairly good agreement. Indeed, to get a more meaningful 
comparison for the self-lensing contribution we normal- 
ize the values for the mass of the luminous components 
to those of the MEGA models (e.g. for their models A, 
M b = 4.4 x 10 10 M Q and M d = 5.5 x 10 10 M©) and use 
a more broadened disk {H = 1 kpc). In Tab. we report 
the obtained results for our models (reference and boxy, 
now labelled A) with the same bulge and disk mass as in 
MEGA models A, for two values of the disk scale height 
H = 0.3 kpc and H = 1 kpc. From Tab. El we can see 
that our estimate for the (total) number of the self-le nsing 
events are in agreement with the lde Jong et all l|2006f) pre- 



5 Note that we are considering a total extinction in the r- 
band of 0.36 mag (0.68 mag) for the bulge (disk), irrespective 
of the line of sight. 



Table 8. Distribution of the number of self-lensing events 
with the distance from the M31 center for several models. 
In the last column, the same quantity is given for dark- 
lensing assuming the reference model, /i = 0.5 and /j, = 
0.2. 



Events inside 










the 8 MEGA 


ref. 


box. 


Wee. 


ref. 


fields 










d(arcmin) 


self 


self 


self 


dark 


2 - 5 


3.81 


4.55 


3.92 


0.93 


5 - 10 


2.80 


3.22 


3.43 


1.97 


10 - 15 


1.32 


1.06 


1.73 


2.46 


15 - 20 


0.66 


0.27 


0.88 


2.20 


20 - 25 


0.42 


0.17 


0.55 


1.78 


25 - 30 


0.22 


0.11 


0.28 


1.23 


30 - 35 


0.09 


0.06 


0.12 


0.78 


35 - 40 


0.04 


0.03 


0.05 


0.33 



diction only when considering more extreme (maximal) 
parameters for the disk component 6 . 

Nevertheless, at variance with Ide Jong et alJ l|2006l) 
we do not conclude that all the 14 events detected by 
the MEGA collaboration can be explained by self-lensing 
only. Indeed, the spatial distribution of the events occur- 
ring inside the 8 MEGA fields, given in Tab.|S]for several 
models (reference, boxy and WeCapp) and shown (nor- 
malized to unity) in Fig. [5] for both self-lensing (refer- 
ence and boxy) and total (self+dark) lensing (reference), 
clearly shows that the distribution with the distance from 
the M31 center of the self-lensing events hardly can be 
reconciled with the MEGA data. Indeed, as seen in Fig. 
an excess of events with respect to expectations from 
self-lensing remains at large distance. This conclusion is 
enhanced assuming the boxy model for the M31 bulge. 

A better agreement with MEGA data can be obtained 
if one considers also a dark-lensing (with fj,^ = 0.5 and 
/h = 0.2) contribution. The compatibility between the ob- 
served MEGA event distribution as a function of distance 
from M31 center and the expected one has been evalu- 
ated for both self-lensing and self+dark lens ing hypothe- 
ses 7 . By using the Kolmogorov-Smirnov test l|Press et all 
Il98fj|) we find a K-S probability ~ 0.51 for self+dark lens- 
ing and ~ 0.18 for self-lensing only, thus implying that 
a dark matter contribution to microlensing seems to be 
favored. 



6 As concerns our estimate in Tab. of dark-lensing events 
due to M31 halo, we obtain a larger number of events with 
respect to MEGA expectations (~ 9 events instead of ~ 6 
events for mh = 0.5 and /h = 0.2). However, to describe the 
M31 dark matter halo we are adopting a different density law 
(an isothermal profile truncated at R — 150 kpc), which is in 
any case consistent with the full M31 rotation curve. 

7 The comparison has been done excluding one event from 
the MEGA candidate list (in the exterior region) since we ex- 
pect that at least one of them is due to the contamination of 
background supernovae (see next for more details.) 
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Fig. 5. For the reference (dotted line) and boxy (dashed 
line) models, the (normalized) distribution of the expected 
number of self-lensing events within the 8 MEGA fields 
is given as a function of the the distance from the M31 
center. The same quantity is shown for self+dark lensing 
(thin solid line) assuming the reference model, /h = 0.2 
and /^h = 0.5. For comparison the (normalized) distribu- 
tion of the 14 observed MEGA events is also given (thick 
solid line). 
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Fig. 6. The expected event number within the 8 MEGA 
fields is given as a function of ti/ 2 , for both self-lensing 
(dotted line) and dark-lensing (dot-dashed line) in the case 
of the reference model. For comparison the distribution 
with tx/2 of the 14 observed MEGA events is also given. 
Here we take fh = 0.2 and ^ = 0.5. 



However, we caution that the candidate microlensing 
events could be contamined by variable stars. In partic- 
ular, the events labelled 13 and 14, located in a region 
where the microlensing rate is negligible, might be con- 
taminated by background supernova e(SN) . Indeed, by as- 
suming standard SN rate IjCoxL 120001) and integrating over 



the volume within z rnax ~ 0.4 (the maximum distance at 
which the SN signal-to-noise ratio is at least 3 a above 
the typical baseline of 22 mag arcsec -2 ) we expect about 
one detectable SN in the outer M31 regions during the 
observational MEGA campaign. 

The distribution of the expected number of events with 
the time scale ti/ 2 is shown in Fig. H3 for the reference 
model and /Zh = 0.5. From this figure, one can see that 
self-lensing and dark-lensing events have almost the same 
tx/2 distribution. Therefore, the t x / 2 event distribution is 
not particularly useful to discriminate the nature of the 
14 MEGA events, at least for a MACHO mass value near 
0.5 M q-, (see also discussion on this point in llngrosso et alJ 
(2006)). The excess of long duration events in the MEGA 
data suggests also a contamination by other variable ob- 
jects. 

We emphasize that our analysis shows that hardly all 
14 MEGA events can be due to self-lensing events by M31 
stars. On the other hand, given the few events detected 
up to now, it seems also premature to derive an estimate 
of the halo dark matter fraction in form of MACHOs. 
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